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1. Introduction 



The singular behaviour of QCD amplitudes is an important ingredient in understanding the perturba- 
tive structure of quantum field theories. In general, when one or more final state particles are either 
soft or collinear, the amplitudes factorise into the product of an amplitude depending on the remaining 
hard partons in the process (including any hard partons constructed from an ensemble of unresolved 
partons) and a factor that contains all of the singularities due to the unresolved particles. One of the 
best known examples of this type of factorisation is the limit of tree amplitudes when two particles 
are collinear. This factorisation is universal and can be generalised to any number of loops [2]. 

Such factorisation properties play a dual role in developing higher order perturbative predictions for 
observable quantities. On the one hand, a detailed knowledge of the structure of unresolved emission 
enables phase space integrations to be organised such that the infrared singularities due to soft or 
collinear emission can be analytically extracted [3, 4, 5]. On the other, the collinear limit plays a 
crucial role in the unitarity-based method for loop calculations [6, 7, 8, 9]. 

For next-to-leading order (NLO) predictions one of the key ingedients is knowledge of the the single 
unresolved configurations where a gluon is soft or where two particles are collinear. The universal be- 
haviour of tree amplitudes in the single soft and collinear limits is well known and has been extensively 
discussed in the literature (see for example Refs. [10, 11, 12, 13, 14]). Antenna factors that describe 
both limits simultaneously have also been derived [15] and subsequently employed in the construction 
of parton level NLO predictions for e + e~ — > 4 jets [16]. At next-to-next-to-leading order (NNLO) one 
encounters, for example, tree amplitudes with double unresolved configurations where two gluons are 
soft [17, 18] or when three particles are collinear [19, 20, 21] as well as the single unresolved limits of 
one-loop amplitudes [22, 6, 23, 24, 25, 26]. Currently, there is significant effort to extract the infrared 
singularities from such processes [27, 28, 29, 30, 31, 32, 33, 34, 35, 36] and to combine them with the re- 
cently computed two-loop amplitudes for parton-parton scattering [37, 38, 39, 40, 41, 42, 43, 44, 45, 46], 
massless Bhabha scattering [47], light-by- light scattering [48, 49] as well as the gluonic production of 
photons [50] and 7* — > 3 partons [51, 52, 53] to provide NNLO estimates of the dominant QCD 
processes in electron-positron annihilation and hadron-hadron collisions. 

One of the recent highlights in perturbative QCD is the succesful computation of the three-loop 
contributions to the Altarelli-Parisi kernels that determine the scale evolution of parton densities. 
These splitting functions have recently been computed in a very impressive series of papers by Moch, 
Vermaseren and Vogt [54, 55, 56] using traditional methods and superseding previous approximations 
based on a limited number of moments [57, 58, 59]. However, an alternative method using the collinear 
factorisation properties of QCD amplitudes has been proposed by Kosower and Uwer [60]. At this 
order one encounters contributions from tree amplitudes with four collinear particles [61], one-loop 
graphs with three collinear particles [62] and two-loop graphs with two collinear particle [1]. 

The two-loop splitting functions describing the time-like splitting of one particle into two collinear 
massless particles are the subject of this paper. They have been studied by Bern, Dixon and Kosower [1] 
who used the unitarity method of sewing together tree (and one- loop) amplitudes to construct two-loop 
amplitudes in the limit where two gluons are collinear both in the case of QCD and N = 1 and N = 4 
super Yang-Mills. In this approach many of the problems associated with using the light-cone gauge 
in traditional Feynman diagram computations of the splitting functions are avoided. In adition, they 
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have checked that the infrared singularities of the splitting function agree with the general formula of 
Catani [63, 64] and provided an ansatz for the form of the non-trivial colour structures that appear 
at 0(1/ e) therein. 

In this paper, we exploit the universal factorisation properties of QCD to extract the two-loop splitting 
functions for all parton combinations by taking the collinear limits of existing two-loop calculations 
of four-point helicity amplitudes [65, 66, 67, 68, 69] and relating them to the two-loop amplitudes for 
vertex graphs. For example, the quark-gluon splitting function is obtained through to 0(e ) by taking 
the collinear limit of the helicity amplitudes for 7* — ► qqg computed in the 't Hooft-Veltman scheme [70] 
and relating it to the two-loop amplitudes for 7* — > qq. Similarly, two-loop helicity amplitudes 
for H — > ggg and H — ► qqg in th elarge m$ limit [71] give access to the gluon-gluon and quark- 
antiquark splitting functions. The master integrals for these processes involve planar and non-planar 
two-loop boxes that have been written in terms of two-dimensional harmonic poly logarithms [72, 73] 
(2dHPL) using differential equations [74]. In the collinear limit, these 2dHPL collapse to harmonic 
polylogarithms of a single variable [75] that can be re-expressed as Nielsen poly logarithms. 

Our paper is organised as follows. In Section 2, we define some basic notation for the three- and 
four-particle processes. Although the full amplitudes factorise, it is often convenient to work at the 
level of colour ordered amplitudes [76, 77, 69, 78, 79, 80]. For the processes at hand, the colour 
structure is particularly simple and the colour ordered amplitudes are straightforwardly obtained by 
stripping away the colour matrices. The collinear singularities of the full amplitude are then obtained 
by restoring these colour matrices. The collinear limit of the colour ordered amplitudes is described in 
Section 3 where we define the splitting functions in terms of the momenta and helicities of the collinear 
particles and also specify how the collinear limit is approached. The procedure for renormalising the 
splitting functions is also described. We turn to the quark-gluon collinear limit in Section 4 and treat 
the gluon-gluon and quark-antiquark collinear limits in Sections 5 and 6. In each case we show how the 
general helicity structure of the four-particle amplitude behaves in the singular limit and give results 
for the tree, one-loop and two-loop splitting functions. In the case of gluon-gluon splitting, we recover 
the QCD results of Ref. [1] in the 't Hooft-Veltman scheme where external particles are treated in 
4-dimensions while the particles internal to the loop are in (4 — 2e)-dimensions. We examine the limits 
as w — > and w — ► 1 in Sec. 7 and show that the soft behaviour of the two-loop splitting functions 
is well behaved and consistant with the expected behaviour. As a by product of our calculation, we 
provide an expression for the two-loop soft splitting function in Section 8. Finally, our main results are 
summarized in Section 9 and accompanied by concluding remarks. The one-loop singularity operators, 
I\ in Catani's language, for each of the processes studied here are detailed in the appendix. 

2. Notation 

We consider the soft and collinear limits of two-loop matrix elements by taking the limits of known 
matrix elements. In general, we consider the decay of a massive object X into three partons, 

X(q) — ► a(p a ) + b( Pb ) + c( Pc ), (2.1) 

where X is a virtual photon or Higgs boson and a, b and c are massless partons. Specifically, we have 

7*^q + Q + 9, (2.2) 
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H^g + g + g, (2.3) 
H^q + g + q. (2.4) 

In the collinear limits, the matrix elements factorise onto the two particle final states, 

X(q)^A(p A )+B(p B ), (2.5) 

or equivalently 

7*^Q + Q, (2.6) 

H — >G + G. (2.7) 

The kinematics of these processes is fully described by the invariants 

Sab = (Pa + Pbf , S ac = (p a + p c f , S bc = (p b + p c f , (2.8) 

and 

Sabc = (Pa + Pb + Pcf = (PA + PB) 2 = SAB, (2.9) 

which fulfil 

Sab + S ac + Sbc = S a bc = SAB- ( 2 - 10 ) 

The unrenormalised amplitude for process V is a vector in colour space and can be written in terms 
of a colour-stripped amplitude, 

\M V ) = C V A V . (2.11) 
For the processes at hand the overall colour structures are given by, 

(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 

where upper and lower indices represent colour indices in the adjoint and fundamental representations 
respectively. Note that the colour ordering of the two- and three-particle amplitudes is trivial here. 
The normalisation Tr(T A T B ) = 8 AB /2 is used throughout. While the explicit calculations presented 
here will be derived from the processes of Eqs. (2.2) -(2.7), the universal factorisation properties of 
QCD amplitudes ensure that the singular limits will be applicable to processes with more particles 
and more complicated colour structures. 

The colour ordered amplitudes have a perturbative expansion in terms of the colour stripped i-loop 
amplitude \A^), 







Sab, 




■*qqg = 


rriC 

1 abi 


C H - 


+GG = 


Tr(T A T B 


C H - 


^999 = 


j?abc 


C H - 


^qqg = 


rpC 

1 abi 



Arp = C-p 



(2.17) 



- 3 - 




a, X, 



A, A., A 



Figure 1: The collinear behaviour of a tree-level amplitude. Particles a and b form (a slightly offshcll) particle A. 
The singular behaviour is encapsulated in the right hand vertex and is denoted by Split A ^ ab '^ (—Xa, a Xa ,b Xb ). 



where C-p contains the overall couplings such that 



a 



V 47ra\/4vra s , 



■y'^-qqg 
Ch~*GG = Clj 

Cn^qqg = Ci\/4ira s , 
and where C\ is the effective Higgs-gluon-gluon coupling. 



(2.18) 
(2.19) 
(2.20) 
(2.21) 
(2.22) 



3. The collinear limit 



When two particles become collinear, colour ordered amplitudes factorise in a universal way [6, 23, 25]. 
If we focus on particles a and b with helicities A a and becoming collinear so that s a t — > and forming 
particle A with helicity A^, then at tree-level we find that, 

A^(a x \b x \...) ^ Yl Split A ^ ab ^\-X A ,a x \b Xb )A^(A XA ,...). (3.1) 

Aa=± 

The factorisation of the (n+ l)-particle colour ordered tree amplitude into a n-particle tree amplitude 
multiplied by the splitting function is illustrated in Fig. 1. The singular factors are all contained in 
Split A ~* a6 '(°)(— Aa, a Xa ,b Xb ) which behaves as l/^/s^. The momentum of the slightly off-shell particle 
A links the two factors, 

PA = Pa + Pb- 

while the helicity A^ is summed over. However, for some helicity combinations, such as the case of a 
positive helicity gluon splitting into two negative helicity gluons, the tree splitting function vanishes. 

Similarly, at one- and two-loops, the colour stripped amplitudes factorise as follows, 

A^(a x «,b x \...) ^ SpKt A - ab ^(-\ A ,a x «,b x »)A(°HA XA ,...) 

\ A =± 

+ Y SpUt A ^ abm {-\A,a Xa ,b Xb )A {1 \A XA ,...), (3.2) 

X A =± 

A^(a x «,b x »,...) ^ Y Split A ^ ab ^(-X A ,a x ^b x »)A^(A XA ,...) 

Aa=± 
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Figure 2: The collincar behaviour of a one-loop amplitude. 



+ Split A ^ ab '^(-X A ,a Xa ,b Xb )A W (A XA ,...) 



x A =± 



+ Split A ^ ab ^ \-XA,a Xa ,b Xb )A {2) (A XA ,...), 

Aa=± 



(3.3) 



where SplhV 4 ~~ >ab,( '*- ) (— A^, a Xa , b Xb ) denotes the i-loop splitting function. These factorisations are illus- 
trated in Figs. 2 and 3. 

In addition to the explicit helicity and particle-type dependence, all of the splitting functions also 
depend on the momentum fraction w carried by each particle, 

k 2 



Pa 



Pb 



(1 - w)p A + kj_- 



n 



1 — w 2pA ■ n ' 



wpa - k± - 
k 2 , 



A:1 



n 



Sab = ~- 



w 2pA • n 
-0, 



(3.4) 



w(l — w) 

where k± is the transverse momentum and n is an arbitrary lightlike momentum such that 2Q ■ k± = 
2n • k± = 0, and implicitly on the dimensional regularisation parameter e. 

The tree and one-loop splitting functions have been known for some time [13, 6, 23, 25], however the 
form of the two- loop splitting function is only known in the case where a gluon split into gluons [1] . 

When the tree splitting function does not vanish, it is customary to factor it out, 

Split A ^ ab ' (1 \-X A ,a x %b Xb ) = r^~ >a6 ' (1) (-A J 4,a Aa ,5 Ab )Split^ a6 '(°)(-A J 4,a Aa ,6 Ai '), (3.5) 

such that all of the infrared singularities lie in the ratio factor rg^ ab '^ 1 \ This factor contains all of 
the non-trivial w dependence as well as the scale violating factor, 

'47r^ 2s 

Sab 
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Figure 3: The collinear behaviour of a two-loop amplitude. 



that modifies the V-y/s^T, behaviour of the tree-level splitting function. r A ^ ab '^ a i so contains the 
infrared poles produced by soft and collinear particles circulating in the loop. Explicit formulae, valid 
to all orders in the dimensional regularisation parameter are given in Refs. [23, 25]. 

Similarly, at two-loops, when the tree-level splitting function does not vanish, the singularities are 



factored into r 



A^ab,(2) 



Split A ^ a6 <( 2 )(-A A , a Xa ,b Xb ) = rg^ ab ' (2) (-X A ,a Xa ,b Xb )Sp\it A ^ ab '^(-X A ,a Xa ,b x >'). (3.6) 

Catani [63] has shown how to organize the infrared pole structure of the two-loop contributions renor- 
malized in the MS scheme in terms of the tree and renormalized one- loop amplitudes. Motivated 
by the structure of Catani's formula and by the results for planar amplitudes in maximally super- 
symmetric Yang-Mills theory [81], we write the singular behaviour of the unrenormalised r A ^ ab '^ as 
(dropping the helicity labels), 



„i4-a6,(2), x _ 



W - ~ 2 (r. 



A-t-o6,(l), 



(0 + 



-«7 



ere 



cr(2e) 



+rp ab ' {2) ' fin (e) + 0(e), 



(3.7) 



where the Euler constant7 = 0.5772 . . . and 



cr(e) 



r(i-e) 2 r(i + e) 

r(l - 2e) 



(3.8) 
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The first two terms contain all infrared singularities from 0(l/e ) through to 0(l/e 2 ) while the 
remaining single poles are absorbed in AH A ~^ ab '^ (e) , 

AH^\e) = {^£j e ( W (l - «,))"* (h?> + ff< 2 > - Hf - (3 K + ft) . (3.9) 

For splitting functions involving only two collinear particles, there are no colour correlations in 
AH A ^ ab >( 2 \e). Here, H q (H q ) and H g are the usual constants appearing in Catani's renormalised 
formula, 

+ (^-54 + Y^J T R N F C F , (3.10) 

^ 2) = |^iVl + T R C F N F - g + § ) TrN f Ca 
Cs , 5 , 11 _ 2 \ ^ 2 



while the constants (3q, Pi and .K" are 

_ 

o 

fl 17C1 - 10CUr fl JV F - 6C f TrN f 

Pi = g , {6.16) 

For SU(N) gauge theory, CU = iV, C F = (iV 2 - l)/2/N and Tr = 1/2. Note that while the pole 
contribution to AH A ~* ab 'W = (hP + - - f3 K + ft) /e is fixed, there is considerable free- 
dom in defining the 0(e) contribution to the factor multiplying the pole. . Different choices affect the 
finite remainder. We choose this particular factor since it regulates the s a b — ► 0, w — > and w — > 1 
limits of the splitting function. 

The remaining finite contributions are contained in rg~~* ab '^'^ in (e) . As we will demonstrate, explicit 
calculation shows that Eq. (3.7) fixes the infrared pole structure of the splitting function exactly, 
leaving only r ^^ ab, ( 2 )'^ m ( e ) to be determined. 

3.1 Ultraviolet renormalization 

The splitting functions have the perturbative expansion, 



Split^ 6 = v 7 ! 



ira. 



SpUt ^a6,(0) + ( ^\ gpli^-a^l) + ( p\ S plit^ a6 '( 2 ) + 0(a 3 s 



(3.15) 



The renormalization of the splitting function is carried out by replacing the bare coupling a s with the 
renormalized coupling a s (/j, 2 ), evaluated at the renormalization scale fi 2 



a s nfS e = a s (/i 2 )/x 2e 



(3.16) 
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where 

S e = (4vrre-% 

and /j,q is the mass parameter introduced in dimensional regularization [82, 83, 70] to maintain a 
dimensionless coupling in the bare QCD Lagrangian density. 

Applying the renormalisation procedure, we have 

Split A-,ab,(0),R = Split A^at,(0) j 

Split A-,ab,(l),R = g-1 Split A^a6,(l) _ Po Split A-a6,(0) ? 

Split ^a6,(2),R = s -2 Split A^a6,(2) _ ^-1 Split ^a 6 ,(l) _ ffh _ ||\ Split ^a6,(0) _ (317) 

Equivalently, for the amplitudes where the tree splitting function does not vanish, 

A^ab,(l),R _ Q -i A^ab,(l) Po 

r s - ^ r s ~ Ye 



„A-a6,(2),fl _ s -2 r A^ab,(2) _ Wo g-i A->ofe,(l) _ / A _ 3/3| . )x , 
e 5 2e e 5 I 4e 8e 2 ' 



Applying this transformation to the unrenormalised splitting ratios, we find that, 



r s 



A^ab,(2),R {e) = 1 ^a 6 ,(l),* (e) N 2 _ ^ r A-H*,(l)^ (c) + e-^(e) ^ + ^ r ^a6,(l), R(2e) 



+Aff ^«6,(2),R( e ) + r A^ab,(2),R,fin^ + Q ^ ^ ^ 

where, 

± H A^ab, { 2),R {e) = £TZ^£) ( f^T) ^ (u;( l _ „,))"* (#(2) + ^(2) _ ( g ^ 

together with a modification of the finite remainder, >ab '( 2 )>- R '-^ n ( e ). xjp i Q trivial gamma function 
coefficients, which only modify the finite remainder, Eq. (3.19) has the same structure as that given for 
the renormalised gluon splitting function in Ref. [1] and also Catani's formula for the general structure 
of renormalised two-loop amplitudes [63]. Throughout the remainder of the paper we will work with 
unrenormalised amplitudes. 

4. The quark-gluon splitting amplitudes 

The unrenormalized amplitude \M) for the virtual photon initiated processes, (2.6) and (2.2), can be 
written as 

\M^ Q ) = S AB V< i %{Q;Q), 

\M r ^ m ) =Tl b V^S,(q;g;q), (4.1) 

where represents the lepton current and S M denotes the colour ordered hadron current. Using the 
Weyl-van der Waerden spinor notation (see Ref. [69]), the hadronic current S M is related to the fixed 
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helicity currents, S Ab , by 



J AB\ 



L 



cj^ b S Ab (Q-;Q+) 



S M (Q+;Q-) 
S M (Q-;Q+) 
S M (g+;</A;g-) = Rj i: f 2 V2a^ B S AB (q+; gX;q-), 
%(q-,g\;q+) = Ll f2 V2a£ B S AB (q-,g\;q+). 



(4.2) 
(4.3) 
(4.4) 
(4.5) 



Labelling the particle momenta by their type, and dropping terms that vanish when contracted with 
the lepton current, we find that 



S Ab (Q+;Q-) = AQ a Q b 



c / , _ \ , s (Iad1 D 1b Rf , 9 AD q D QB , , Qcb9°9A 
S AB (q+;g+;q-) = a{y,z) y _ w ^ + 0{y,z) tttvtt^ + 7U/, *0 



(qg)(gq) 



(qg)(gq) 



(qg)(gq)* ' 



(qg)*(gq)* 



(qg)*(gq)* 



(qg)*(gq) 



(4.6) 
(4.7) 

(4.8) 



In these expressions, (ij) = (i — \j+) and [ij] = (i + |j— ) where |iA) is the Weyl spinor for a massless 
particle with momentum i and helicity A. The spinor products are antisymmetric and satisfy [ij] = 
-(ij)* and (ij)[ji] = 

The currents with the quark helicities flipped follow from parity conservation: 



S AB (Q-,Q+) = (S BA (Q+;Q-)y, 
S AB (q-;gX;q+) = (S BA (q+;g(-X);q-))*. 



(4.9) 



The coefficients A and Q (Q = a, (3, 7) contain all of the integrals over loop momenta and ultimately 
determine the infrared structure of the amplitude. A depends only on the invariant mass of the 
virtual photon sab, while Q depends on y = s qg /s qgq and z = s gq /s qgq . Each coefficient is a linear 
combination of master loop integrals [72, 73] with coefficients that are rational functions of the scale 
invariants and the spacetime dimension. Expansions of the renormalised Q around e = are given in 
Ref. [51, 52, 53]. 

Taking the quark-gluon collinear limit corresponds to y = s qg /s qgq — > 0, with the quark and gluon 
carrying momentum fractions (1 — w)Q and wQ respectively. When the quark and gluon both have 
positive helicity we find that, 



S^ B (g+;5+; Q-) 



Q\\9 



1 



w (qg) _ 



<x(y,z) (l-w)+P(y,z) 

y-»0 



w 



QaQ 



A^B- 



(4.10) 



Similarly, when the quark has positive helicity but the gluon has negative helicity, we find, 

1 — w 



S AB (q+;g-;q-) 



q\\g 



Vw [qg] [ 



a(z,y) 



y->0 



Q a Qb- 



(4.11) 



We see that in both cases we recover the two-particle helicity structure of Eq. (4.10). 



-9- 



In the antiquark-gluon collinear limit, z = s g q/s qg q — > and the antiquark and gluon carrying mo- 
mentum fractions (1 — w)Q and wQ respectively. When the antiquark has negative helicity then, 

1 



S AB (q+;g-;q-) 
S AB {q+;g+;q-) 



Q\\9 

\/w 

q\\g 1-W 



w 



z->0 



w (gq) 



a(y,z) 



z->0 



QaQb- 



The helicity amplitude coefficients have the perturbative expansions, 

2 



A - C 7 »_qq 



* m +(s) A(,, + (s) * M +<x» 



and 



for fi = a, /3, 7. 



n (0) (y,^) + (g) n (1) (v,*) + (£)^ (2) (^) +<?( 



Qa<3b> (4.12) 
(4.13) 



(4.14) 
(4.15) 



4.1 The tree- level quark-gluon splitting amplitudes 

At leading order, 



a 



(0) 



(y, z) = (y, z) = 1 and ^4 (0) = 1. 



(4.16) 



Bearing in mind that the splitting function defined in Eq. (3.1) relates the colour ordered amplitudes, 
comparing Eqs. (4.10) and (4.11) with Eq. (4.6), we can immediately read off the tree-level quark-gluon 
splitting functions for colour stripped amplitudes, 

V2 



SplitG-™(°)(-, g +,0+) 
Split^'(°)(-,g+ <T) 



Vw (qg) ' 
-V2(l-w) 



(4.17) 



Vw [qg] 

The factors of y/2 are remnants of the definition of the helicity current in Eq. (4.4). In addition, the 
soft gluon singularity as w — > is explicit. 

Splitting functions where Q has negative helicity are obtained by complex conjugation according to 
Eq. (4.9), 

-y/2 



Spht^(°)(+,g-, 5 -) 
Split Q ^°\+, q-,g + ) 



Vw [qg] ' 
y/2(\-w) 
w (qg) 



(4.18) 



Similarly, the tree splitting functions for an antiquark Q to split into a gluon g and antiquark q are 
given by, 



Split "^ ) (+,(?+,<?-) 
S plit^-,(o) ( _ 



^ , SpUt^<°>(+,g-,n = - ^-y (4.19) 



w 
V2 



w (gq) 



, Split^(°)(-, 5 + ,^ 



Vw (gq) 
y/2(l - w) 
w [gq] 



(4.20) 
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4.2 The one-loop quark-gluon splitting amplitudes 



The unrenormalised one-loop splitting functions denned in Eq. (3.5) are presented in Ref. [25] in a 
variety of schemes to all orders in e, and we reproduce them here in the 't Hooft-Veltman regularisation 
scheme (Sr = 1) for the sake of completeness, 



<9-5S,(l), 



-cr(e) 



T,q ,g 

2\ e 



Q^ot,(i), ± 
' c> vi ? y 



5 



[2e 2 
iV 2 + 1 



m=l 



9 

LL 



— W 



N J 4(1 - 2e) 



-cr(e) 



47T/X 2 



iV 
2? 



oo 

Li m 



?Tt=l 



1 — W 



As usual, the polylogarithms Li n (w) are defined by 

f w dt 

Li n (w) = / — Li n _i(i) 
Jo f 

Li 2 (w) = - / -ln(l-t). 
Jo t 



iV2 Lim 



iV2 Lim 



for n = 3, 4, . . . 



1 — w 



(4.21) 



(4.22) 



(4.23) 



We have checked that our procedure of taking the collinear limit of the explicit one-loop amplitudes 
written in terms of 2dHPL reproduces these expressions through to 0(e 2 ). Note that, unlike Ref. [25], 
in our notation the gluon always carries a momentum fraction w. Note also that because of our different 
normalisations, the one-loop splitting function differs by a factor of N/2 compared to Ref. [25]. 



4.3 The two-loop quark-gluon splitting amplitudes 

Explicit calculations for the unrenormalised two-loop coefficients, that independently satisfy the pole 
structure predicted by Catani (see Appendix A) can be found in Ref. [52]. Using FORM [84] 
and MAPLE to take the limits of the 2dHPL, we find that the unrenormalised splitting functions 
rg^"'^ 2 \e) have a pole structure determined by Eq. (3.7). The finite remainders are given by, 
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5. The gluon-gluon splitting amplitudes 

The gluon-gluon splitting amplitudes can be extracted from the Higgs decay into gluons in the limit 
where the top quark is very heavy and an effective Higgs-gluon-gluon vertex is induced. Labelling the 
gluon momenta alphabetically we write the unrenormalised amplitude in terms of the single colour 
ordered amplitude, 

\M h ^ab) = SabA h (A,B), (5.1) 

where, 

A H (A+,B+) = A H [AB] 2 , (5.2) 
A H (A + ,B-) = 0. (5.3) 

The amplitudes where the first gluon has negative helicity can be obtained by parity. 
Similarly, we write the Higgs to three gluon amplitude as, 

\M H ^abc) = f abc A H (a,b,c). (5.4) 

Again, there is only one colour ordering and the two independent helicity amplitudes for the three 
gluon decay are [85, 71], 

A H (a+, 6+ c+) = -a H (y, z) ±= f« bc (5.5) 

V2 {ab){bc){ac) 

1 \ac\ A 

A H (a+,b-,c+) = p B (y,z) -= - (5.6) 

V 2 [ab\ [bc\ [ca\ 

where y = s ab /s abc and z = s ac /s abc . 

The other helicity amplitudes are obtained by the usual parity and charge conjugation relations. For 
the present purposes, the only useful one is, 

A H (a-,b-,c + ) = p H (z,y) -L ff 4 , (5.7) 

V2 {ab} {be) (ca) 
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The limit where gluons a and b are collinear corresponds to y = s a b/ s afec — * 0- We assign a momentum 
fraction (1 — w) to particle a. In this limit, we see that, 



A H (a + ,b + ,c + ) 
A H (a + ,b-,c+) 
A H (a~ ,b~ ,c + ) 



allfe 1 



\[2 y / w(l-w)(ab) 

a\\b 1 (1-w) 2 



®H(y,z] 



y-,0 



[AB] 2 , 



- — ; nl^M )[AB] 2 , 



allb 1 



<a6) 



^^(l-t«)[a6pV4 



-Ph(z,v) 



[ABf 



Note that the terms in round brackets are evaluated in the collinear limit, y — > 0. 
The two-particle helicity amplitude coefficient A# has the perturbative expansion, 



^ff = Ch^gg 



4? + (£K*+S) 2 ^+^ 



>2tt/ " V2vr. 

Similarly, the three-particle helicity amplitude coefficients an, and (5h are given by, 

2 



ftff (2/, z ) = C H ^ggg 

for n# = a H ,PH- 



ng) (y> z) + (g) n«(y, z) + (g) ng>(y, z) + O(o») 



5.1 The tree- level gluon-gluon splitting amplitudes 

At leading order 

«2W) = /4°W) = 1 and AP = 1, 
and we can immediately read off the tree- level gluon-gluon splitting functions, 

Split G ^< (0 V,& + ) ' 



and by symmetry, 



SpUt?-™ (0) (a+6-) = - 
Split G ^ s ' {0) (a-,6-) = 0, 



S pIitf- +flfl>(0) (o-,6 + ) 



y/2 y/w(l - w) {ab) 
1 (1-w) 2 



V2 y/ w {\ - w) [ab] ' 



w 



V2 yX 1 - w) [ab] ' 

The amplitudes where a gluon with negative helicity splits are obtained by parity, 



Split^ 99 ' (0) (a\& A >) = (Split G ^' (0) (a-\r A <>))\ 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



and recalling that [ab]* = —(ab). 
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5.2 The one-loop gluon-gluon splitting amplitudes 

As in the quark-gluon case, the one loop splitting functions are presented in Ref. [25], and we reproduce 
them here in the 't Hooft-Veltman regularisation scheme (Sr = 1) for the sake of completeness. For 
the amplitude that vanishes at tree-level we have, 



Spht_ yyA '(a ,b ) = cr 



2\ « 



w(l — w) 



-s ab J (l-2e)(2-2e)(3-2e) 
1 (ab) 



(N -Ne- N F ) 



V2 y/w{l-w)[abX' 



(5.15) 



while the other splitting functions are, 



tt;(l — it;) 



(l-2e)(2-2e)(3-2e) 



(N - Ne- N F ) 



(5.16) 



G^gg,(l) 



{T,a T ,c 



cr(e) 



4vr^ 2 V N 



SabJ 2e2 



\ / m=1 \ / 



(5.17) 



Note that Eqs. (5.16) and (5.17) are symmetric under w <-> (1 — it;). Note also that because of 
our different normalisations, the one-loop splitting function differs by a factor of N/2 compared to 
Ref. [25]. We have checked that we reproduce these expressions through to 0{e 2 ). 

5.3 The two-loop gluon-gluon splitting amplitudes 

Using explicit calculations for the unrenormalised two-loop coefficients for the H — > gg and H — > ggg 
decays, we find that the unrenormalised splitting function r^ 9g ^ 2 \e) for the helicity combinations 
that are non- vanishing at tree level have a pole structure determined by Eq. (3.7) with finite remainders 
given by, 
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G^>gg,(2),fin, + 



(-,</ + ,<T) = 



,r2 / 3 , / x 3 ! \ 11 , 5 . 523 \ 

+/?o ^-Li 3 (io) + i Li 2 (iy) ln(u>) + J ln(u;) + i ln(l— w) (2 + j ln(l-u/) 

55 25 65 \ 
+ 24 C3 + 48 C2 "2T6j 

"^ (T^) ( Li sW " Cs - ^ In(w) (Li 2 H - C2) 
-AT(7V - 7V F ) f ^J^ (li 3 M - Ca - \ M«0 (Li 2 («;) - (2! 



-Li 2 (l-w) + 2C2 + 



1 / ln(w) 



6(1 -u>) 2 I v ; * / 12 \(l-w) 
Finally, for the amplitude that vanishes at tree- level we find, 
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(5.20) 



We have checked that the expressions presented here agree with those given in Ref. [1]. 
6. The quark-antiquark splitting amplitudes 

The quark-antiquark splitting amplitudes can be extracted from the Higgs decay into quarks and 
gluons in the large mt limit, which factorises onto the H — > GG amplitudes in the collinear limit. The 
single colour ordered amplitude is obtained from the full amplitude by, 



\M H ^ q qg) =T c ab A H (q,q,g). 



(6.1) 



There is one independent helicity amplitude for the H — > qgq decay which is given by [85, 71], 

1 \qq? 



Ah(q ' ,q ,g ) = ih(v,z) 



V2 



(6.2) 
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where y = s qq /s qqg and z = s qg /s qqg . All other helicity amplitudes are either zero due to helicity 
conservation or can be related by charge conjugation and parity. 

In the collinear quark and antiquark limit, y = s qq /s qqg — > 0. When the quark carries momentum 
fraction (1 — w)) we see that, 



A I + — +\ 1 W I 



3/— 



[AB]< 



Note that the terms in round brackets are evaluated in the collinear limit, y — ► 0. 

The three-particle helicity amplitude coefficient 7#, has a perturbative expansion given by, 



m{y,z) = c H - 



•qqg 



(6.3) 



(6.4) 



6.1 The tree- level quark-antiquark splitting amplitudes 

At leading order 

(o) n 



(6.5) 



Together with the tree amplitude for H — > GG in Eq. (5.2), we can immediately read off the tree-level 
quark-antiquark splitting function, 



SpUt^ M "' (0) (g+g-) = 



1 w 



V2 [qq] ' 

6.2 The one- loop quark-antiquark splitting amplitudes 



(6.6) 



As in the quark-gluon and gluon-gluon case, the one loop splitting functions are presented in Ref. [25]. 
In the 't Hooft-Veltman regularisation scheme (5r = 1), 
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(6.7) 



Note that because of our different normalisations, the one- loop splitting function differs by a factor of 
N/2 compared to Ref. [25]. We have checked that we agree with this expression through to 0(e 2 ). 



6.3 The two-loop quark-antiquark splitting amplitudes 

Using explicit calculations for the unrenormalised two-loop coefficients for the H — > gg and H — > 
qqg processes, we find that the unrenormalised splitting function rg~^ q ^ 2 \e) has a pole structure 



- 18 - 



determined by Eq. (3.7) with a finite remainder given by, 
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7. Limiting Properties 



It is useful consider various limiting cases of the splitting functions in order to get a more physical 
picture of what is going on. 
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Table 1: Leading behaviour of splitting functions 



Process 


Split* ) 


Split* 1 ) 


Split* 2 ) 


G + - g-g- 





v 7 ! — tc 


v 7 ! — w 



Table 2: Leading behaviour of splitting functions 
as w — > f 



In Tables 1 and 2 we show the leading behaviour of the splitting functions as a polynomial in w for 
io — ► and (1 — w) for w — > 1. aj,6j,Cj and are the i-loop coefficients multiplying the leading 
term. Both limits correspond to the production of a soft particle. When the particle is a gluon, and 
the helicity of the hard particle is inherited from the parent, the tree splitting function is singular 
(Q + — ► q + g ± , G + — > as w -> and G + — ► as «; — > 1). In each of these cases, the one- and 

two-loop splitting functions produce the universal constants a\ and 02, indicating that the soft limit 
is independent of the particle type and the helicity of the gluon. In the soft quark (antiquark) limit, 
there is no polynomial enhancement. 

Furthermore, when the helicity of the hard particle is flipped (as in the Q + — ► q + g~ and G + — > g + g~ 
splittings asi/)-> I) there is a tree- level suppression factor of at least (1 — w). In the two-loop results 
presented in the previous sections, the finite contributions to the helicity flip processes Q + — > q + g~ , 
G + — ► and G + — > splittings contain apparent quadratic, cubic and quadratic singularities 

as w — > 1. In each case, the singularity is softened by the behaviour of the numerator, such that only 
the G + — > <7 + g~ splitting has a residual singularity of 1/(1 — to). This is of course annihilated by the 
tree-level helicity flip suppression factor of (1 — w) 3 ^ 2 . 

All of these limits show that the soft behaviour of the two-loop splitting functions is well behaved and 
consistent with the expected behaviour. 



8. The two-loop soft splitting function 



In the limit that a gluon becomes soft, the colour ordered amplitudes factorise in a similar way to 
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Eqs. (3.1), (3. 2) and (3.3). At tree level we have the usual result, 

\M^Ha^^,c x %...)) b ^S^\a^,c)\M^ 1 (a^,c x %...)) 
The tree level soft splitting amplitudes are given by, 

V2\ac] 



S(°)(a,6-,c) = — 



(8.1) 



(8.2) 



[ab][bcY 

The soft amplitudes are independent of the helicities and particle types of the neighboring legs, a and 



c. 



The factorisation of one-loop colour ordered amplitudes in the soft limit proceeds as, 

\M^{a x «,b x \c x %...)) ^ s(°)(a,6\c)|A4£ 1 2 1 (a^,c A V..)> 

+5( 1 )(a,6^, C )|Mi°2 1 (a A «,c^,...)). 



(8.3) 



The onedoop soft factor has been found, to all orders in e, to be [25, 26] where it is expressed in terms 
of the treedevel soft factor, 



SW(a,b ± ,c)=rW ft (a,b ± ,c) 5(°)(a, ^.c), 



cr r iv 



(a, 6 ,c) = -- 



soft y 



r(l + e)r(l-e). 



where, 

Twodoop amplitudes have a similar soft limit, 

b^O „„,), , \ , . f (2) , A a A 



(8.4) 

(8.5) 
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(8.6) 



As in the onedoop case, it is convenient to factor out the overall singularity as the tree soft factor, 



S {2 \a,b ± ,c) 



,(2) 
soft 
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We can use the twodoop quark-gluon and gluon-gluon splitting amplitudes presented in Sees. 4 and 

5 to recover the soft limit. Compared to the collinear limit defined by Eq. (3.4), the soft limit is 

obtained by taking w small. In terms of invariants, w — ► Sf, c /s a i). Expanding the explicit expressions 
for r G*^9<y,(2)(_ )ff + )5 +) ? r G-flff,(2)(_ >fll + jflr ) ^ r Q^qg,(2)( T ^ q ±^ g ±j and r Q->qg,(2) g ± ^ gT ) in the 

to — > limit, we find the universal twodoop soft splitting factor to be, 
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Note that although the soft limit has simple QED-like factorisation properties for colour ordered 
amplitudes, the soft limit of the colour dressed amplitudes is rather more complicated. Nevertheless, 
the procedure for colour dressing is straightforward [26, 1]. 

9. Conclusions 

In this paper we have studied the universal behaviour of two-loop amplitudes in the limit where two 
massless external partons become collinear using explicit calculations of 7* — > qqg and H — > 3 partons. 
While we work with two-loop amplitudes in the specific case, our results generalise to processes 
involving more external particles. Each of the amplitudes has a trivial colour structure, so our results 
naturally apply to colour stripped (or ordered) amplitudes. 

We find that the pole structure of the unrenormalised two-loop splitting amplitudes is given by 
Eq. (3.7). After renormalisation, this is in agreement with the pole structure predicted by Catani [63] 
for generic renormalised two-loop amplitudes. A similar structure has been found for the gluon split- 
ting function in Ref. [1]. 

The finite remainders for the two independent quark-gluon splitting functions are given in Eqs. (4.24) 
and (4.25), while similar results for the three different helicities needed to describe gluon-gluon splitting 
are given in Eqs. (5.18), (5.19) and (5.20). In this latter case, we find complete agreement with the 
results obtained in Ref. [1] using a more general approach involving the unitarity method of sewing 
together tree (and one-loop) amplitudes to directly construct the two-loop splitting functions. Finally, 
the finite two-loop remainder for the one independent quark-antiquark splitting function is given in 
Eq. (6.8). 

Our results are valid for the time-like splitting of colour ordered amplitudes in the 't Hooft-Veltman 
dimensional regularisation scheme. The procedures for extending their range of validity into the 
space-like region or dressing with colour are straightforward [1] . 

We have studied the limits where one of the split partons is soft and, although there are apparent 
additional polynomial singularities as w — > 1 in some of the splitting functions, find that the logarithms 
and polylogarithms multiplying the explicit poles conspire to reduce the degree of the singularity by 
two powers. The net result is that the soft splittings which induce a helicity flip on the hard parton 
remain suppressed. 

For the processes where the helicity of the hard parton is conserved, we find that the soft limit of the 
radiated soft gluon is universal and is independent of the helicity or particle type of the hard parton. 
This universal factor is the two-loop soft splitting function and we give an explicit expression for it in 
Eq. (8.8). 

The splitting amplitudes presented here are another tool in the armoury of perturbative QCD. They 
can serve as a check on calculations of two-loop QCD amplitudes with more external legs. Together 
with the triple collinear limits of one-loop integrals [62] and the quadruple collinear limits of tree 
amplitudes [61], they may also be useful in rederiving the recently calculated [55, 56] three-loop 
contributions to the Altarelli-Parisi kernels [60]. Finally, they will provide one of the ingredients 
necessary for NNNLO computations of jet cross sections. 
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A. Infrared factorisation of two-loop amplitudes 



Catani [63] has shown how to organize the infrared pole structure of the two-loop amplitudes renor- 
malized in the MS scheme in terms of the tree and renormalized one-loop amplitudes. In particular, 
the infrared behaviour of the one-loop coefficients is given by 
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(0)\ 
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where |Afp^) is the em renormalised z-loop amplitude. Similarly, the two-loop singularity structure is 
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where the infrared singularity operators 1$ and 1$ are process but not helicity dependent. Each 
of these operators produces terms of 0(e) and the division between finite and non- finite terms is to 
some extent a matter of choice. Adding a term of 0(1) to is automatically compensated for by a 
change in I^ 2 \ 

For the virtual photon initiated processes, (2.6) and (2.2), 
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Note that on expanding Sy, imaginary parts are generated, the sign of which is fixed by the small 
imaginary part +i0 of Sij. 

For the Higgs processes, (2.7), (2.3) and (2.4), we have 
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The two- loop singularity operator Ji, is given by, 




r(l - 2e) f/3o 
r(l-e) V e 




In principle, ii^ contains colour correlations between the final state particles. However, for the 
processes at hand, these correlations vanish and each external coloured leg in the partonic process 
contributes independently to H v 



where n g is the number of external gluons, n q is the number of external quarks and anti-quarks for 
the process V and H q and H g are defined in Eqs. (3.10) and (3.11) respectively. 
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